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Abstract 

Optically thin two-temperature accretion flows may be thermally and viscously 
stable, but acoustically unstable. Here we propose that the 0-mode instability of 
a cooling-dominated optically thin two-temperature inner disk may explain the 23- 
day quasi-periodic oscillation (QPO) period observed in the TeV and X-ray light 
curves of Mkn 501 during its 1997 high state. In our model the relativistic jet 
electrons Compton upscatter the disk soft X-ray photons to TeV energies, so that 
the instability-driven X-ray periodicity will lead to a corresponding quasi-periodicity 
in the TeV light curve and produce correlated variability. We analyse the dependence 
of the instability-driven quasi-periodicity on the mass (M) of the central black hole, 
the accretion rate (M) and the viscous parameter (a) of the inner disk. We show 
that in the case of Mkn 501 the first two parameters are constrained by various 
observational results, so that for the instability occurring within a two-temperature 
disk where a = 0.05 — 1.0, the quasi-period is expected to lie within the range of 8 
to 100 days, as indeed the case. In particular, for the observed 23-day QPO period 
our model implies a viscosity coefficient a < 0.28, a sub-Eddington accretion rate 
M ~ 0.02 MEdd and a transition radius to the outer standard disk of tq ~ 60 r^, 
and predicts a period variation 5P/P ~ 0.23 due to the motion of the instability 
region. 
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1 Introduction 



Strong variability is one of the common observational properties of TeV emit- 
ting Active Galactic Nuclei (AGNs) [10]. In many cases, the highly variable 
high-energy gamma-rays and the X-rays appear to be correlated with no time 
delays evident on day-scales, suggesting that the 7-rays may result from in- 
verse Compton upscattering of lower energy photons. The first TeV blazar, 
for example, that was observed simultaneously in multiple bands from radio 
to TeV gamma-rays is Mrk 421. The first campaign, conducted in 1994 [52] 
on Mrk 421, revealed correlated variability between the keV X-ray and TeV 
gamma-ray emission. The gamma-ray flux varied by an order of magnitude 
on a timescale of 2 days and the X-ray flux was observed to double in 12 hr. 
On the other hand, the high-energy gamma-ray flux observed by EGRET, as 
well as the radio and UV fluxes, showed less variability than the keV or TeV 
bands. Another multiwavelength campaign on Mrk 421 performed in 1995 re- 
vealed another coincident keV/TeV flare [6|63] . The UV and optical bands 
also showed correlation during the flares. With the detection of TeV emission 
from Mrk 501 [S2], several multiband campaigns were organized on Mrk 501 to 
verify whether such a behavior is a general property of TeV-emitting blazars 
or whether it is unique to Mrk 421. 

The gamma-ray blazar Mkn 501, detected as a strong TeV emitter in 1995 
|52] . is one of the closest {z = 0.0337) and brightest BL Lacertae objects. 
Historically (i.e., prior to 1997), its spectral energy distribution (SED) u 
resembles that of X-ray selected BL Lac objects, having a peak in the extreme 
UV-soft X-ray energy band p^fSTj . Earlier optical observations of Mkn 501 
have shown variations of up to l.'"3 and polarized emission up to Popt — 7% 

During its active state in 1997, where Mkn 501 was monitored in the 2-10 keV 
X-ray band by the all sky monitor(ASM) on board RXTE and in the TeV 
energy band by several Cherenkov telescope groups [T|9lll7f32ll51ll56] . both X- 
rays and TeV gamma-rays increased by more than 1 order of magnitude from 
quiescent level PH5] . Analysis of the X-ray and TeV data showed, that the 
variations were strongly correlated between both bands, yet only marginally 
correlated with the optical UV band PI17II32|47] . While the synchrotron emis- 
sion peaked below 0.1 keV in the quiescent state, in 1997 it peaked at 100 keV. 
This is the largest shift ever observed for a blazar |38]. Earlier investigations 
of the 1997 April flare in Mkn 501 [IH], showed that its origin may be related 
to a variation of 'jmax together with an increased luminosity and a flattening 
of the injected electron distribution. 

During the 1997 high state, the X-ray and TeV light curves displayed a quasi- 
periodic signature [2^301151] . A detailed periodicity analysis, based on the 
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TeV measurements obtained with all Cherenkov Telescopes of the HEGRA- 
Collaboration and the X-ray data with RXTE, was performed by Kranich et 
al. [5UIIM] and more recently also by Osone [IS] including Utah TA data. The 
results indeed strongly suggest the existence of a 23 day periodicity, with 
a combined probability of p ~ 3 x 10~^ in both the TeV and the X-ray 
light curves covering the same observational period [30|3T] . Rieger and 
Mannheim (2000) have shown that the origin of these QPOs may be related 
to the presence of a binary black hole in the center of Mkn 501. While this may 
well be possible, we explore here an alternative scenario, where the observed 
QPOs are related to accretion disk instabilities. 

In a seminal paper, Shapiro et al. (1976) (SLE) have pointed out that a 
hot (Compton cooling-dominated) optically thin two-temperature accretion 
disk may be present in the inner region of the standard optically thick disk 
|59] . whenever the SLE inequality a^^/^M^M*"'^^^ > 0.6 is satisfied, where 
= M/(10^^gs-^), = M/(3M0), M is the accretion rate of the disk, 
M the mass of the central black hole, and a the viscosity parameter, con- 
strained by the model to lie within the range 0.05 < a < 1. Later work P9]|50] 
has shown that the SLE configuration might be thermally unstable (although 
less than the standard disk) if the simple standard viscosity prescription is 
employed. However, relatively small changes in the viscosity law can already 
ensure a stable configuration [19], in particular, if stabilizing effects of a disk 
wind are fully taken into account. A kind of SLE two-temperature disk struc- 
ture may thus well exist in the inner region of BL Lac type objects [66ll67f8] and 
provide a possible explanation for the X-ray and TeV variability phenomenon 
in Mkn 501. For, firstly, Compton processes in a inner two-temperature disk 
with electron temperature Tg of about 10^ K (and ion temperature one or 
two orders of magnitude higher) can produce (steep) X-ray power-law spec- 
tra, in contrast to the standard disk model that can only produce emission 
up to the optical-UV band [SSf^ . Secondly, analysis of the linear stability of 
an optically thin two-temperature disk around a compact object shows that 
the disk is subject to a radial pulsational instability (inertial-acoustic mode 
instability) [69]. This possible mode of pulsational over stability, in which ra- 
dial disk oscillations with local Keplerian frequencies become unstable against 
axisymmetric perturbations, occurs if the viscosity coefficient increases suffi- 
ciently upon compression In this case, thermal energy generation due to 
viscous dissipation increases during compression, leading to amplification of 
the oscillations analogous to the role played by nuclear energy generation in 
stellar pulsations. As we demonstrate below, the occurrence of such a type of 
disk oscillation may well account for some quasi-periodic time variability in 
AGNs in a way similar to Galactic black hole candidates [Sll23|40|70j . 



2 Note that no QPOs have been seen by MAGIC during 1998-2000, when the 
source was not very active [2j, suggesting that the QPOs only occur during a very 
active stage. 
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2 Model and Results 



2.1 Model 

Due to the Lightman & Eardley secular instability [3^11^7] of the standard 
geometrically thin, optically thick disk model [SH], a two-temperature disk 
configuration is likely to be present in the inner region of the standard disk. 
As shown by Shapiro et al. [58] , the outer boundary of such a two-temperature 
disk is determined by 

riy\~'{ro) ~ 3.85 x lO^a'/'Mj'^'Ml. (1) 

where tq = rg*^^ is the radius of the outer boundary, Mj is the mass of the 
central black hole in units of 10"^ Mq, M24 is the accretion rate in units of 10^^ 
g s~^ = 0.016Mq/jt, a is the viscosity parameter constrained to be within the 
range of 0.05 to 1.0, and ( = 1 — (^)~^^^ expresses the boundary condition 
that the viscous stress must vanish at the inner edge of the disk. 

The total luminosity of a two-temperature disk, integrated from the inner edge 
to To, is of order 

LrOA4 = ^[l + 2(^ro.)-'/'-18/ro,]M24 (2) 
4 

where ^^0,44 is in units of 10^^ erg s~^. Combining relations (1) and (2), one 
obtains 

-1 

= 2.62 X 10^L^o,44A (3) 

where A = {aM^'^Y^^ . Obviously, since A is relatively insensitive to a, the 
outer radius tq* and the accretion rate M can be uniquely determined, once 
the central black mass M and the observed radiative output ^^0,44 are known. 

Following Shapiro et al. [58] (see also [66]) we assume, that a hot (cooling- 
dominated) optically thin two-temperature region forms within the inner re- 
gion of the standard disk. Hence, the classical outer and middle regions of the 
standard Shakura &: Sunyaev (SS) disk model [SH] describe the outer portions 
of our disk model. The inner region of the cool SS disk, whose outer boundary 
lies at the point where the radiation pressure Pr matches the gas pressure 
Pg, extends inwards to the radius tq, where Pr = 3Pg- The radius marks 
the outer boundary of the two-temperature inner region, which then extends 
inward to the innermost radius 3rg for a non-rotating black hole. In short. 



l + 2(^)-3/^-i^ 

6 rn*. 



^21/16 



^^0* \-l/2 

^ 6 
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the assumed disk configuration is assumed to be identical to the standard SS 
model for r > Tq, but described by the two-temperature structure equations 
for r < tq. The disk thickness, the density, and the ratio ^ for the equilibrium, 
two-temperature inner disk then satisfy the following equations: 

h = 5.29 X lO^i A/^/^'M^'f ^ 

p= 1.14 X 10-"M7-'/'M2-//'«3/4^-1/4^-9/8 ^ ^^-3 (5) 

h/r = 0.316 X M~'/'' M%'' a-'^'' r^^s . (6) 

Typically, the emission spectrum of a hot (Tg ~ 10^ K) inner disk is brem- 
strahlung unless a soft photon source is present, in which case the spectrum 
may become a power-law. In our hybrid model a fraction of the optical-UV soft 
photons from the outer SS disk is assumed to be intercepted by the hot SLE 
inner disk and Comptonized to form the X-ray power-law spectrum, similar 
as in [581166] . If Comptonization is not saturated, as possible in the presence of 
a copious soft photon source (such as a cold SS component), the soft photons 
will be upscattered into a power-law distribution over some energy range. The 
resulting disk spectrum above the input soft energy Eg = hug is approximately 
of the form 

oc 1^-°^ exp(-/iz//A;Te) , (7) 

where = (2.25 + 4/[y (1 + /t)])^/^ _ 3/2 > 0, /t = 2.5 6 + 1.875 e\l - 6), 
9 = /cTg/meC^ and y ~ 1 for unsaturated Compton cooling, i.e., the spectrum 
resembles a power-law with index ax for E > Es and E ^ kTg, and shows 
an exponential cut-off for E ~ kTe [58|66] . The corresponding SLE disk SED 
uFi, peaks at Ep = {1 — ax) kTe, i-e., usually at around some tens of keV, 
with Ep being higher for smaller indices. In particular, for kT^ = 200 keV 
(Te = 2.3 X 10^ K) one may have ax - (2.25 + 1.85 /y) - 1.5 G [0.4, 0.8] for 
y E [0.6,1.3]. 

As first shown by Kato (1978), accretion disks can undergo (radial-azimuthal) 
pulsational instabilities and thereby cause quasi-periodic variability. If a fluid 
element in the disk is perturbed from its equilibrium position, it oscillates 
in the radial direction with the epicyclic frequency, as the difference between 
gravitational and inertial (Coriolis and centrifugal) force appears as a restoring 
force. Due to this property, axisymmetric perturbations can propagate in the 
disk as inertial acoustic waves. As shown by Kato, a viscous disk can be pul- 
sationally unstable against this type of oscillations under certain conditions: 
By studying the stability properties of an optically thin disk to axisymmet- 
ric, local (wavelength A -C radial size of disk) and nearly radial {vr ^ Vz) 
oscillations, he found that the disk becomes unstable if the coefficient of the 
viscosity increases sufficiently rapidly with density and temperature, and that 
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the frequency of oscillations essentially corresponds to the angular frequency 
Q of the disk, independent of wavelength. More recently, Wu (1997) has stud- 
ied in detail the radial pulsation instability of optically thin two-temperature 
accretion disks, demonstrating that in a geometrically thin, cooling-dominated 
two-temperature disk the acoustic 0-mode (outward-propagating) is always 
unstable. It has been argued, that pulsational overstability might explain some 
of the QPO phenomena observed in cataclysmic variables (CVs) and Galactic 
black hole systems [5|ll|40ll43j . and be important for a proper understand- 
ing of periodic variability in AGNs [23]. As we show in the the next section, 
Mkn 501 may indeed represent a promising AGN candidate, where pulsational 
overstability becomes apparent. 

Numerical simulations of pulsational instability of geometrically thin, optically 
thick (one-temperature) a disks show that inertial acoustic waves are excited 
and propagate both inwards and outwards periodically, immediately growing 
to shock waves and resulting in a time- varying local accretion rate exceeding 
the imposed input value M by one or two orders of magnitude . Global 

nonlinear calculations of such systems indeed demonstrate that (i) the net 
mass flow changes sign by a significant amount according to the direction of the 
oscillatory flow, although the disk seems to maintain its ability to transport on 
average the initially imposed M [27], and that (ii) the oscillations also induce 
small relative changes in other variables (e.g., temperature or surface density) 
[l3]. In particular, the luminosity variations caused by overstable oscillations 
in such systems are expected to be (only) at the percent level, even if local 
variables such a mass accretion rate and radial velocity change significantly 
[TT|27|43j . Unfortunately, to our knowledge, no simulation of pulsational over- 
stability for a SLE type, optically thin, cooling-dominated two-temperature 
disk configuration has been performed up to now. This makes it difficult to 
draw solid (quantitative) conclusions for the SLE the simulations per- 

formed so far consider physical environments quite different from that of a 
hot, optically thin, gas-pressure supported two-temperature accretion disk. 
Global simulations of optically thin, two-temperature disk configurations are 
certainly required to remedy this problem. Nevertheless, it seems very likely 
that significant changes in local accretion rate during an oscillatory period 
represent a genuine qualitative feature associated with pulsational overstabil- 
ity. If so, then substantial quasi-periodic changes in the total disk luminosity 
might be expected for the SLE case (cf. [SH], eqs.[l] and [21]). 

In principle, the observable period of the pulsational instability in a viscous 
disk is of order [7D] 



= 5.78 X 10"^ (1 + z) M7r2/2a-i[days] , (8) 
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where = ^Jr'^/GM is the local Keplerian timescale, r,,, in units of 
is the radius at which the pulsation instability occurs, and z is the redshift 
of the source. Obviously, since the disk structure depends significantly on the 
radius, the period of the pulsational instability is different from one radius to 
another. Yet, by using eq.(3) we can calculate the radius ro, at which the disk 
becomes unstable to radial pulsation for a given luminosity, and so derive an 
upper limit on the range of possible periods. In fact, these upper limits may 
offer a very useful guide, since numerical disk simulations suggest that the 
oscillations may always be trapped near the boundaries of the discs [27^6] . 

2.2 Application to Mkn 501 

The 7-ray blazar Mkn 501 was the second source after Mkn 421 to be detected 
at TeV energies in 1995 by the Whipple observatory [52]. In 1997 the source 
came into much attention when it went into a remarkable state of strong 
and continuous flaring activity, becoming the brightest source in the sky at 
TeV energies [T|2|9|T7|22f51j . BeppoSAX observations during a strong out- 
burst in April 1997 showed, that the spectrum was exceptionally hard {a < 1, 
Fi, oc over the range 0.1-200 keV, indicating that the X-ray power out- 
put peaked at 100 keV or higher energies [H]. These observations implied 
that the peak frequency has moved up by more than two orders of magnitude 
(persistent over a timescale of ~ 10 d) and that the (apparent) bolometric 
luminosity of Mkn 501 has increased by a factor of ~ 20 compared to previ- 
ous epochs pHl32|48|57|64j . Optical observations on the other hand indicated, 
that the source was still relatively normal and thus suggested that the varia- 
tions are confined to energies above ~ 0.1 keV ^7\. Further observations with 
BeppoSAX in April-May 1998 and May 1999 during periods of lower TeV 
flux showed that the peak frequency has decreased to ~ 20 and ~ 0.5 keV, 
respectively [6^ . 

A comprehensive analysis of the temporal characteristics of the gamma-ray 
emission from Mkn 501 in 1997 showed that the TeV emission varied signifi- 
cantly on time scales of between (5-15) hours [1]. If this short time variability 
is via inverse Compton processes related to accretion disk phenomena (see 
§ 3), we can estimate an upper limit for the black hole mass in Mkn 501 by 
taking the variability timescale to be of order of the Keplerian orbital period 
at the innermost stable orbit r^s, suggesting 



for a non-rotating black hole where r^s = 3rg = 2^%-, and thus a black hole 
mass of M < (3.8 — 11.5) x IO'^Mq for the (5-15) hour timescale observed. Note, 




(9) 
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that this mass upper hmit may be somewhat higher, if a pseudo-Newtonian 
potential is employed, in which case one may find M < (5.7 — 17.2) x 10'^ Mq. 
Eq. ([H]) essentially assumes that the observed short-time variability is domi- 
nated by processes occurring close to the innermost stable orbit. This seems 
justified as long as the associated timescale {t ~ 3.6 x 10^ sec in the lab. 
frame) is larger than (i) the transverse light crossing time of the source as 
measured in the lab. frame, i.e., U ~ r6/[rbc] < t, where ~ z/Vb with 
Ffe ~ 15 the bulk Lorentz factor and z ~ 5 x lO^r^ the distance from the 
central engine, cf. [12], and (ii) the characteristic shock acceleration timescale 
^acc ~ 10'^ rgyr/(rfec) < ti [nS], where is the electron gyroradius for intrinsic 
magnetic field strengths ~ 0.01 G and comoving electron Lorentz factors of 
~ 10^ (see below), which seems indeed to be the case. Based on the above 
noted, as well as related considerations in [51], we adopt a black hole mass 
of M = 9 X 10'' Mq, for which Vg ~ 2.66 x 10^^ cm, as fiducial value in our 
calculations below. 

Detailed observations of Mkn 501 during the 1997 active phase showed that 
its X-ray emission was highly variable, although the soft X-ray fiux (up to a 
few keV) did not change dramatically [52|34|41|H5] . BeppoSAX observations 
in 1997 give (2-10) keV fiuxes in the range (1.69 — 5.24) x 10~^° erg cm~^ 
s~^, with an average value of 2.50 x 10"^'^ erg cm~^ s~^ [H], corresponding 
to a mean (bolometric) luminosity of < Lx >= 6.2 x 10^^ erg s~^, assuming 
Ho = 72 km/s/Mpc and a fiat Universe with Qm = 0.27. If the observed 
periodicity in the (2-10) keV band is caused by a pulsational disk instability, 
a non-negligible fraction / < 1 of this X-ray luminosity has to be produced 
by the inner disk, the other part being provided by the relativistic jet. Each 
surface of the disk then produces a characteristic (2-10) keV luminosity of 
-f'2-iokev = 3.1 X 10^^ / erg s~^. Simulations of pulsational instability discussed 
at the end of § 2.1 indicate that the oscillations are accompanied by periodic, 
large amplitude variations in local accretion rate (exceeding the average M by 
over an order of magnitude) and possibly trapped near the outer edge of the 
disk. Simple quantitative modelling then suggests that the total luminosity of 
a SLE type disk may vary by up to a factor of a some few, i.e., / may be 
as small as ~ 0.1. Using the observed soft X-ray power-law index — 0.7 
[53|41f47j . we may estimate the observationally required, total integrated SLE 
disk luminosity L^o by calculating the ratio r] = i^2-iokcv/-^ro of the soft X-ray 
to total disk luminosity from the emitted disk spectrum (see eq. [7]), which 
for kTe ^ 200 keV gives rj ~ 1/5. This suggests that the total (average) SLE 
disk luminosity is of order L^o ^ 16 x 10^^ / erg s~^. 

As shown below, low / values are more preferable in the hybrid SLE model and 
we will thus henceforth adopt a total SLE disk luminosity < Lro >= 1.75 x 10^^ 
erg s~^ (/ = 0.11) for explicit calculation (see however also Fig. [Tjfor the more 
general case). Given the constraints on the viscous parameter a and using 
M = 9.0 X IO^Mq, we may then obtain the following results for the transition 
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radius iq*, the maximum pulsational period P and the ratio -, namely, for 
a = 0.05 we have: ro* = 97, P = 100 days, ^ = 0.55; while for a = 1.0 
one finds: ro* = 130, P=8 days, ^ = 0.09, respectively. In both cases the 

accretion rate is similar, i.e., M24 ^ 2.7, corresponding to M ~ 0.02MEdd, and 
consistent with the condition (see § 1) for a SLE configuraton to exist. We 
note that for BL Lac objects the transition radius between an inner optically 
thin two-temperature and an outer standard disk has been estimated to lie 
within the range tq* — (40 — 150) for M / M-^dd = 0.01 [8], which seems well 
consistent with the values derived above. Based on these results, we can also 
roughly estimate the luminosity Lss ~ associated with the SS disk 

component at r > Tq: Employing the thin cold disk temperature relation 
T ~ 6.3 X lO5mi/^(lO^M0/M)^/4(rg/r)3/4 - I.3 x 10^ K with r = tq ~ 50rg 
one finds that Lgs ^ 1.8 x 10^^ erg/s, and that the emission is maximized at 
frequencies around /~3/cT//i^8x 10^^ Hz. These findings imply a Compton 
energy enhancement factor A considerably larger than Lro/Lgs — 10. 

Spectroscopic data for Mkn 501 suggest, that the normal bolometric (photo- 
ionization-equivalent) luminosity of the disk component driving the NLR 
photoionization is of order a few times 10^^ erg/s, corresponding to a photon 
number per second that can ionize hydrogen of Nii{H) ~ Lf,/(/ii^i) ~ 10^^ s~^, 
where ui = 13.6 eV/h. Here, Lf, has been estimated based on the Lh-L{Ha) 
relation, i.e., logL^ = 1.176 \ogL{Ha) —4.91 [53], with an observed (mean) 
narrow Ha line luminosity for Mkn 501 of ~ 10^^ erg/s [61], which gives L5 ~ 
2 X 10^^ erg/s (with substantial scatter). We note that this value is very close 
to the one estimated by Barth et al. (2002) using emission line measurements 
for four nearby PR I radio galaxies. Given the properties of the cold SS disk 
component derived above, the main part of the ionizing photons per second is 
expected to come from the SLE component (unless its existence is only coupled 
to active source stages). The number of photons emitted by the SLE inner disk 
per second, which can ionize hydrogen is given by N{H) = L^/[hi']di', 
where L^, is the specific SLE luminosity (cf. eq. [7]). Using the values adopted 
above, i.e., ~ 0.7, kT^ ^ 200 keV and La-iokeV = 5.8 x lO^V/ erg s'^ 
(p ~ 0.3 and p = 1 for the quiescent and active source stage, respectively, cf. 
[iT]). we can find N{H) ~ 4.5 x lO^^p/ s"\ which for p ~ 0.2 and / ~ 0.1 is 
comparable to Nb{H), and may thus qualify the choice of / above. 

As mentioned in § 1, detailed periodicity analysis of the X-ray and TeV light- 
curves of Mkn 501 during its high state in 1997 has provided strong evidence 
for a 23-day period in both energy bands [5Ull31|45j . As shown above, such a 
period falls within the range of possible pulsational periods, suggesting that 
it may be well explained by the pulsational instability of a two-temperature 
SLE disk. In particular, by using P = 23 day and eqs. (3) and (7), we can 
place an upper limit on the allowed viscosity parameter and the expected 
transition radius. Por the characteristic values derived above, one obtains a < 
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4e+44 4.5e+44 



Fig. 1. Outer SLE radius tq*, accretion rate M^ot = -^/AfEdd and upper limit for 
the viscosity coefficient a < Oc as a function of the luminosity of a two-temperature 
SLE disk. Parameters M = 9 x 10^ Mq and Pobs = 23 d have been employed in the 
calculations. 

0.28, corresponding to iq* = 115 and ^ ~ 0.20. The more general case is 
shown in Fig. ([1]), where the dependency of these results on the total SLE 
disk luminosity has been calculated. 



The changing SLE disk luminosity will lead to a variable seed photon field for 
Compton upscattering to the high energy regime by the nonthermal particles 
in the relativistic blazar jet. The disk radiation entering a moving plasma blob 
from behind appears de-boosted by a factor ~ F^ in the comoving frame of the 
outflowing plasma. It will nevertheless dominate over possible (quasi-isotropic) 
disk photons that are scattered by diffusive gas or clouds into the direction of 
the source, provided the distances z to the central engine is smaller than 



^ 7 X 10^ 




0.001\i/2 
) ' 



g 1 



(10) 



where -Rsc is the distance and t^c the mean scattering depth of the electron 
scattering cloud [IS] . This usally places the emission site at distances between 
hundred (to avoid 77 absorption) and some thousand r^. 

As shown by Dermer and Schlickeiser (1993), scattering of X-ray (~ 0.1 keV) 
photons entering from behind will generally take place in both the Thomson 
and the Klein-Nishina regime for a large range of (lab. frame) soft photon 
energies eph, i.e., for energies in the range TbmeC^/'j < eph < iTb'jmeC^, where 
Ffe ~ 15 is the jet bulk Lorentz factor. An full treatment of the scattering pro- 
cess is thus a complicated matter and beyond the scope of the paper. However, 
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we can still get some order-of-magnitude insight into the associated photon 
luminosity by approximating the scattering process by its Thomson limit: The 
maximum value of the scattered photon energy e^^ in the lab. frame is given by 
^ 0.5 7^ eph [I3]. Upscattering of soft X-ray photons (~ 0.2 keV) to the TeV 
regime thus requires comoving electron Lorentz factors of 7 ~ 10^. Scattering 
will always take place in the Thomson regime for (lab. frame) soft photon en- 
ergies eph < rfemeC^/7 ~ 0.08 keV [2], suggesting that our approximation is 
indeed not too bad if we take an analogous reduction of the Thomson cross- 
section by a factor of two or so into account. For photons entering from behind 
the single inverse Compton power in the comoving blob frame for scattering 
by a relativistic electron with Lorentz factor 7 then becomes, cf. |15j . 

Pco.p~-aT^^. (11) 



Suppose, that electrons are continuously accelerated at a shock front into a 
characteristic non-thermal (isotropic) differential power law particle distribu- 
tion in the comoving frame of 72(7) = no7~^ for 10^ ~ 7i < 7 < 7max with 
7max ~ 2 X 10^ roughly given by the balance between acceleration and external 
Compton cooling for B ~ 0.05 G and z ~ 500 r^. Due to 7max oc L^q' , an 
increase of the SLE disk field will lead to a decrease in the maximum Lorentz 
factor. As particle are advected downstream and no longer efficiently accel- 
erated, the distribution steepens to n{'y) oc 7""^ for momenta where particles 
have had sufficiently time to cool, resulting into a broken power law electron 
distribution with a break energy of 7;,. We may estimate 7^ by setting the 
(comoving) external Compton cooling time scale equal to the (comoving) dy- 
namical time scale ~ z/{rbc), which gives 7;, ~ 2 x 10^. For typical values 
no7i ~ lO^cm"'^, e.g., [17f33] . the associated observable Compton VHE lumi- 
nosity can thus be estimated from Lobs ~ rft-Pcomp-^eq^^ where iVb ~ ?T-(7b)7b 
and AV ~ An{z/rb)^rb is the volume scale of the emitting blob measured in 
the comoving frame. This gives a Compton luminosity 

^^^^ ^ (50^) (1^) iw^) ' ^^^^ 



which is of the order of the observed VHE luminosity above 1.5 TeV during 
the 1997 high state [T|TF] . This suggests that changes in the disk photon field 
driven by a pulsational instability may indeed be responsible for the observed 
modulation of the TeV emission. During the quiescent state on the other hand, 
the X-ray flux is typically about ten times smaller (cf. §2.2). Hence, even if a 
SLE type disk is present in the quiescent state (see however below), the ex- 
pected TeV contribution due to inverse Compton scattering of disk photon is 
much smaller, perhaps even swamped by synchrotron self-Compton emission. 
The 1997 observations of Mkn 501 showed that the observed TeV spectrum 
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gradually steepened with energy starting at around 3 TeV, approaching a log- 
arithmic spectral {i' FJ) slope above 5 TeV of ~ 1.7 ± 0.33 ± 0.6 including 
systematical and statistical errors PES] . However, as shown for example by 
Konopelko et al. (2003), this curvature at around 3 TeV is hkely to be caused 
by gamma-ray absorption in the intergalactic infrared background field, sug- 
gesting that the peak in the intrinsic (de-absorped) spectrum due to increasing 
Klein-Nishina effects is around 8 TeV. In principle, it appears possible that 
due to the strong angle dependence of the scattered flux on the angle of the 
impinging soft photons, photons from the outer SS disk part, i.e. from disk 
radii ~ z, could make an important inverse Compton contribution as well. 
Redoing the calculations presented in [15] for lab. frame angles of the incident 
photons of ~ 7r/4 gives a total photon energy density of the soft photons in the 
comoving frame of ~ (1 — /3b/-\/2)^ Upj^ where Wp^ g ^ Lss(r0/z)/(47r2;^c) 
is the corresponding energy density in the lab. frame. As the SLE contribu- 
tion Mph^, entering from behind, appears de-boosted by a factor of ~ 4r^ 
[T5] compared to its lab. frame value, the ratio of comoving energy densities 
becomes 



using Fb = 15, L^o/Lss ~ 10 and z/vq ~ 10, indicating that the SLE part 
dominates during the high state even if we assume that the SS disk is not 
truncated at radii of ~ 500 Vg. Note however, that the latter may well be the 
case if a binary black hole system exists in the center of Mkn 501 as argued 
for by several authors (see discussion). 



3 Discussion 

The X-ray to TeV spectra of TeV blazars have been often interpreted within a 
one-zone SSC model, e.g. [^16ll42f62|[M] . However, if a two-temperature disk 
structure is present in some of these sources, the real situation may be more 
complex as the X-ray radiation from a two-temperature disk, for example, may 
represent an additional, non-negligible source of seed photons for the inverse 
Compton scattering to TeV energies, in particular during active source stages, 
likely to be associated with changes in accretion history. Hence our model as- 
sumes, that a pulsational instability occurring within a two-temperature disk 
leads to observable, periodic variations of its X-ray radiation field. Part of this 
periodically modulated X-ray emission will enter the jet and (in addition to di- 
rect synchrotron photons) serves as seed photons for Compton-upscattering to 
TeV energies similar as in [13111^ . Our model thus takes it that both, a direct 
synchrotron self Compton and an external Compton contribution are relevant 




(13) 
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for modelling the SED of Mkn 501, the seed photons for external Compton- 
upscattering consisting of both, the infrared-optical seed photons from the 
(quasi-steady) SS disk component and the variable X-ray photons from the 
two-temperature disk component. A related, but more simpler scenario, as- 
suming the seed photons for inverse Compton scattering to be provided by a 
(direct) synchrotron plus an quasi-steady flux component comparable to the 
observed infrared-optical flux, has been proposed by Plan et al. (1998) in order 
to account for the different degrees of SED variations of Mkn 501 at X-ray and 
sub-X-ray energies during the April 1997 outburst (see also pTpi] ). Detailed 
analysis indeed suggests that one-component SSC models cannot fit both the 
April 1997 SEDs and the lightcurves from X-ray to TeV, and that (at least) 
an additional, moderately variable low energy component contributing in the 
energy range between 3-25 keV is required [33)1^ . A similar conclusion seems 
to hold for the strong X-ray outburst observed in July 1997 [31]. Interestingly, 
observations in 1998 also provide evidence for an additional component in the 
optical regime, possible associated with the SS disk component in our hybrid 
disk model (see also [25] for an alternative interpretation) : As shown by Mas- 
saro et al. [H] optical to X-ray data taken in June 1998 indicate that the 
optical spectrum is steep and does not match the low energy extrapolation of 
the X-ray spectrum, hence suggesting the presence of different emission com- 
ponents in the optical and in the X-ray regime as naturally expected in our 
model. 

Monitors in both the X-rays and the TeV emission show evidence for a 23- 
day periodicity during the 1997 high state. As demonstrated above the 23-day 
period in the X-ray light curve may be caused by a pulsational instability in 
two-temperature accretion disk, and via the inverse Compton process result in 
the same periodicity in the TeV light-curve. A pulsational instability occurring 
in a two-temperature disk with transition radius tq ~ (48 — 65) Vg will result 
in a recurrence timescale of 8 to 100 days. 

Based on the observed TeV variability we have employed a characteristic black 
hole mass of 9 x 10'' in our calculations. We note that quite different central 
mass estimates for Mkn 501 have been claimed in the literature, ranging from 
several times 10^ (mainly based on high energy emission properties) up to 
10^ Mq (based on host galaxy observations), see e.g. [71ll2lll8fl9IIM] . However, 
as shown by Rieger & Mannheim [51] uncertainties associated with host galaxy 
observations may easily lead to an overestimate of the central black hole mass 
in Mkn 501 by a factor of three and thus reduce the implied central mass to 
^ (2 — 3) X 10® Mq, a value in fact recently confirmed by an independent 
analysis of central mass constraints derived from host galaxy observations |68j . 
Moreover, as argued by the same authors some of the apparent disagreement 
in central mass estimates may possibly be resolved if a binary black hole 
system exists in the center of Mkn 501, see also [53|54ll65j . similar as in the 
case of OJ 287 [BO]- For example, if Mkn 501 harbours a binary system with a 
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more massive primary black hole of ~ IO^Mq and a less massive (jet-emitting) 
secondary black hole of ~ 10^ Mq, the mass ratio p = m/M would be of order 
0.1, which may compare well with the result p < 0.25 estimated for OJ 287 
[38] . While our characteristic black hole mass employed falls well within the 
above noted range, we note that the SLE pulsational instability model may 
still work successfully, if a higher black hole mass is used. For example, if one 
adopts M = 3 X IO^Mq and P = 23 days, one obtains a < 0.07, tq* = 46 and 
M ~ 0.008MEdd. 

Our analysis is based on a specific disk model (SLE) which is open to ques- 
tions, in particular with respect to its possible stability properties. We note 
however, that a relatively small change in the usually employed viscosity de- 
scription may already lead to a thermally stable configuration [19]. On the 
other hand, it may as well be possible that the SLE configuration represents a 
quasi-transient phenomenon associated with those changes in accretion history 
that probably initiate the high states. An alternative (inner) disk configu- 
ration of interest may be represented by an optically thin two-temperature 
ADAF solution [Him] . Such a configuration can exist for accretion rates 
m = M/M-Edd below a critical rate merit — 0.3 ^ 0.019, where the canonical 
ADAF value of a = 0.25 has been employed [H]. ADAFs are generally less 
luminous than a standard disk, with the typical ADAF luminosity given by 
La — 0.02 (m/a^) M [7T]. Using the constraints above, the possible ADAF 
luminosity for Mkn 501 becomes < 1.3 x 10^'^ erg/s, which is already about 
an order of magnitude smaller than required by the X-ray analysis. This sug- 
gests that - at least during its high state - an optically thin ADAF is not a 
viable option for Mkn 501. 

Based on Eq. (7) we can estimate the variation rate of the period due to the 
motion of the instability 



' 2 rdt 



(14) 



Using Vr 



dr 

at' 



M = AirphrVr and Eqs. (4) and (5), one finds 



6P/P = 0.44 M^r' Mf r:'/\-'/' (15) 

For a = 0.28 the relevant parameters result in 6P/P = 0.23. From the period 
analysis performed by Kranich et al. [30] a 3a deviation in period corresponds 
to 6.67 days. If we take this deviation as the intrinsic variation on the peri- 
odicity (P), then a 6P/P of ^ = 0.29 can be estimated from the results by 
Kranich et al. [30]. As this estimate assumes that the deviation of the period 
is only affected by the motion of the instability, while it may in fact be caused 
by more than one effect, our theoretical SP/P should not be greater than the 
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observational results. We conclude that the observed 23-day QPOs in Mkn 501 
might be caused by the instability of a two-temperature accretion disk. The 
model presented here may thus offers an alternative explanation to the binary- 
driven helical jet model of Rieger & Mannheim (2000). Comprehensive compu- 
tational modelling of the pulsational instability in a two-temperature, cooling 
dominated disk will be essential to verify this in more detail. 

Our model predicts that a period correlation in the X-ray and 7-ray should 
always be present during an active source stage, while the period of the QPOs 
may vary as the instability region could change from one high state to the 
other. 
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